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ABSTRACT 

Deep and high resolution radio observations of the Hubble Deep Field and flanking 
fields have shown the presence of two distant edge-darkened FR I radio galaxies, allow- 
ing for the first time an estimate of their high redshift space density. If it is assumed 
that the space density of FRI radio galaxies at z>l is similar to that found in the local 
universe, then the chance of finding two FRI radio galaxies at these high radio powers 
in such a small area of sky is < 1%. This suggests that these objects were significantly 
more abundant at z>l than at present, effectively ruling out the possibility that FRI 
radio sources undergo no cosmological evolution. We suggest that FRI and FRII radio 
galaxies should not be treated as intrinsically distinct classes of objects, but that the 
cosmological evolution is simply a function of radio power with FRI and FRII ra- 
dio galaxies of similar radio powers undergoing similar cosmological evolutions. Since 
low power radio galaxies have mainly FRI morphologies and high power radio galax- 
ies have mainly FRII morphologies, this results in a generally stronger cosmological 
evolution for the FRIIs than the FRIs. We believe that additional support from the 
y/Vmaii test for evolving and non-evolving population of FRIIs and FRIs respectively 
is irrelevant, since this test is sensitive over very different redshift ranges for the two 
classes. 



1 INTRODUCTION 



Radio source counts and the cosmological evolution 
of radio galaxies 

Already in the early days of radio astronomy it was re- 
alised that number counts of radio sources contain impor- 
tant information about the distribution of radio sources as 
a function of redshift. Early results showed that they were 
inconsistent with the predictions of a Steady-state cosmol- 
ogy (Ryle & Clarke 1961), but that they could be explained 
by evolutionary cosmological models in which there were 
many more powerful radio sources at epochs earlier than the 
present. This could be directly inferred from the slope of the 
LogA'^— Logs' curve above Jy levels being steeper than the 
— I power law which the static, uniform, Eucleidian Universe 
must show. At sub-Jy levels the slope changes gradually to 
a law flatter than the — | law as it must to escape the radio 
equivalent of Giber's Paradox. Longair (1966) showed that 
the strong evolution must have been confined to the most 
powerful radio sources, because the observed range in power 
for radio sources was more than 5 orders of magnitude, while 
the turnover in the radio source counts covered less than 2 
orders of magnitude. 

The initial steep slope and its turnover were subse- 
quently accounted for by two distinct populations of radio 
sources. One population consisted of high radio luminosity 
sources undergoing strong cosmological evolution (e.g. with 



comoving number densities about 2 — 3 orders of magni- 
tude higher at large redshifts), the other population was of 
low luminosity sources undergoing no or only little evolution 
with cosmological epoch. Wall (1980) suggested that these 
two populations of non-evolving and strongly evolving radio 
sources corresponded to the Fanaroff & Riley (1974) Class 
I and II galaxies respectively. This idea has been further 
developed by Jackson and Wall (1999), who use the source 
counts at several frequencies and redshift information of the 
3C sample, to link the FRI and FRII radio sources as the 
unbeamed parent populations of BL Lac objects and flat 
spectrum quasars respectively. Willott et al. (2000) also de- 
fined a dual-population scheme, but slightly different from 
that of Jackson & Wall (1999). In their scheme, the low lu- 
minosity population is associated with radio galaxies with 
weak emission lines (both FRIs and FRIIs), and the high 
luminosity population is associated with radio galaxies and 
quasars with strong emission lines (almost all FRIIs). In 
contrast, Dunlop & Peacock (1990) do not make a distinc- 
tion between FRI and FRIIs, in their analysis of the statis- 
tics of steep spectrum radio sources, but allow for a lumi- 
nosity dependent cosmological evolution. 



Fanaroff &: Riley Class I and II radio sources 

Fanaroff & Riley (1974) showed that there is a distinct mor- 
phological difference between radio sources of high and low 
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luminosity. The more powerful sources have their regions 
of highest surface brightness at the ends of a double lobed 
structure (FR II — edge brightened), while the lower power 
objects show a variety of forms in which the highest bright- 
ness occurs near their centres, excluding their cores (FR I 
— edge-darkened). Fanaroff & Riley (1974) found a sharp 
division between the two classes at Pitsmhz ~ 2 x 
Hz-^sr-^ 

Early optical work indicated that there is also a differ- 
ence between the host galaxies of FR I and FR II type radio 
galaxies. Owen and Laing (1989) found that FRII sources 
reside in normal giant elliptical galaxies with absolute mag- 
nitudes near Af, of the Schechter luminosity function, con- 
siderably fainter than first-ranked galaxies in rich clusters, 
while FRI sources reside in galaxies which can generally be 
described as D or cD galaxies. However, more recently it 
has been shown by Ledlow and Owen (1996) that this result 
was caused by the combination of a sample selection effect 
and a strong positive correlation between the FRI/II radio 
luminosity cutoff and the absolute magnitude of the host- 
galaxy, with an increase of the transition luminosity from 
Li.4GHz«102'' W/Hz at Mfl=-21, to Li.4GHz~102'^ W/Hz 
at Mi{=— 24. Due to the small range of radio-luminosity in 
the Owen and Laing sample, this correlation resulted in the 
FRIIs generally residing in lower luminosity galaxies than 
the FRIs. The Ledlow & Owen result suggests that both 
FRI and FRII radio sources live in similar environments, 
but that the properties of the host galajcies may strongly in- 
fluence the morphological appearance of the radio sources, 
at least for those near the division of the two classes. 

Zirbel & Baum (1995) found that for the same to- 
tal radio power, FRII galaxies produce 5—30 times more 
emission- line luminosity than FRI galaxies. In this light, 
Baum, Zirbel & O'Dea (1995) put forth the possibility that 
the FR dichotomy is due to qualitative differences in the 
structural properties of the central engines in these two types 
of sources, like the accretion rate and/or the spin of the cen- 
tral black hole. In contrast, Gopal-Krishna & Wiita (2000) 
point out a class of double radio sources in which the two 
lobes exhibit clearly different FR morphologies. Although 
these objects are rare, their existence supports explanations 
for the FR dychotomy based upon jet interaction with the 
external medium, and appears quite difficult to reconcile 
with the class of explanations that posit fundamental differ- 
ences in the central engine. 

Radio Observations of the Hubble Deep field 

Due to observational limitations, no direct measurements of 
the high-z space density of FRI radio galaxies yet exist. The 
Hubble Deep Field (HDF; Williams et al. 1996) and the sur- 
rounding Hubble Flanking Field (HFF) is the best studied 
area of sky to date. The random region is unbiased, although 
it was chosen for the absence of any bright objects at any 
wavelength. Recently, several groups have observed the HDF 
in the radio regime, using the VLA, MERLIN, WSRT and 
the EVN (see table |l]). This allows a first estimate of the 
high redshift space density of FRI radio sources, which is 
discussed in this paper. 

Richards et al. (1998) have observed the HDF at 8.4 
GHz using the VLA in A, CnB, C, DnC and D configu- 
rations, corresponding to angular resolutions ranging from 



0.3" to 10". They detected 29 radio sources in a complete 
sample within 4.6' of the HDF centre and above a flux den- 
sity limit of 9 /iJy, of which 7 are located in the HDF itself. 
Muxlow et al. (1999, 2001) have observed the HDF at 1.4 
GHz, using the VLA for 42 hours and MERLIN for 18 days. 
They have detected a complete sample of 87 sources in a 
10' X 10' region with flux densities above 40^Jy. These have 
all been imaged with the MERLIN-(-VLA combination to 
produce images with 0.2, 0.3, and 0.5" resolution with an 
rms noise-level of 3.3/iJy. In addition, Garrett et al. (2000) 
have observed the HDF with the WSRT at 1.4 GHz, re- 
sulting in an angular resolution of 15" and a noise- level of 
8 /iJy. They detect 85 regions of radio emission (> 5cr) in 
a 10' X 10' fleld centered on the HDF, with 4 sources not 
previously detected at 1.4 GHz. Garrett et al. (2001) have 
taken deep European VLBI Network (EVN) observations of 
the HDF, with a resolution of 25 milliarcseconds. This has 
resulted in the detection of 3 sources above a flux density 
level of 210 ^Jy- 

The deep optical/radio observations of the HDF and 
HFF indicate that ~ 60% of the faint sub-mjy and /ijy radio 
sources are identifled with star forming galaxies at moderate 
redshifts [z ~ 0.2 — 1), often showing morphologically pe- 
culiar, interacting/merging galaxies with blue colours and 
Hll-like emission spectra (Richards et al. 1998). Another 
20% of the sources are as yet unidentified in the optical 
down to 1=25.5 in the HFFs and 1=28.5 in the HDF. These 
faint systems may be distant galaxies obscured by dust. The 
remaining 20% of the faint radio population seem to be iden- 
tified with relatively low luminosity AGN. 

This paper concentrates on the two brightest radio 
sources at 1.4 GHz, both exhibiting a clear FRI radio mor- 
phology. We will discuss the implications of their presence in 
the HDF and HFFs for the cosmological evolution of these 
objects. Throughout this paper we assume a cosmology with 
Ho=50 km sec~^ Mpc~^ and qo=0.5, for consistency with 
previous studies. 



2 FRI RADIO GALAXIES IN THE HDF AND 
HFFS 

HDF123644+621133 

The radio source HDF123644-I-621133 is the brightest ob- 
ject in the HDF at 8.4 GHz. It has an 8.4 GHz flux density 
of 477 AtJy. Muxlow et al. (2001) give a total 1.4 GHz flux 
density of 1.2 mjy using VLA-I-MERLIN data, while Gar- 
rett et al. (2000) report a total flux density of 1.6 mJy using 
the WSRT. This higher flux density is probably a result of 
the lower resolution of the WSRT observations, containing a 
larger proportion of the total source flux. The radio source 
consists of an unresolved core with a relatively flat spec- 
trum (ai.4-8.4GHz = —0.3 ± 0.2, with f{v) oc !/"), which 
is surrounded by steep spectrum ai.4-8.4GHz = —1.2 ± 0.2 
emission oriented north-south and extending about 30" (Fig 
The radio source shows a prototypical FRI morphology. 
It is optically identifled with a bright -Rab = 20.5 elliptical 
galaxy at z=1.013 (Richards et al. 1998). Assuming the flux 
densities and spectral indices of the core and extended struc- 
ture given above, HDF123644-I-621133 has an observed flux 
density at 178 MHz rest wavelength (88 MHz observed) of 
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Telescope 

VLA A,CnB, C, DnC , 
VLA + MERLIN 
WSRT 
EVN 



Frequency Resolution Flux limit Number of Reference 

(GHz) (arcsec) (A'Jy) Sources 

D 8.4 0.3-10 9 29 Richards et al. (1998) 

1.4 0.2 40 87 Muxlow et al. (2001) 

1.4 15 40 85 Garrett et al. (2000) 

1.4 0.025 210 3 Garrett et al. (2001) 



Table 1. Radio observations of the Hubble Deep Field and Flanking Fields. The columns give the telescope used, the observing frequency, 
resolution, flux density limit, number of detected sources, and the reference to the observation. 




Figure 1. The 1.4 GHz image of the FRI radio sources HDF123644+621133 from VLA-data (left) and of HFF123725+621128 from 
VLA+MERLIN data (right), kindly provided by T. Muxlow (Muxlow et al., in preparation 2001). 



SnsMHz = 24 mjy. This corresponds to a total radio power 
of Pi78MHz = 102'' « W Hz-isr-1. 



HFF123725+621128 

The radio source HFF123725+621128 is the brightest radio 
source in the HDF and HFFs at 1.4 GHz, with a total flux 
density of ~ 6 mJy. The high resolution MERLIN+VLA im- 
age, crucial for a morphological classification, shows it to be 
5" in extent with two high brightness regions (peak bright- 
ness of 200 /ijy/beam) close to its centre, making it an FRI 
radio source (Fig. |l]). The two diffuse extended lobes at the 
east end west side are both slightly bent towards the south 
giving the object the appearance of a wide angle tail (WAT), 
which at low redshift are mostly found in galaxy-cluster en- 
vironments. Richards et al. (1998) quote a spectral index for 



this source of a1.4_g.4GHz ~ —1.0 ± 0.1. Its detection in the 
WENSS survey at 325 MHz (Rengelink et al. 1997) of 25 
mJy is consistent with this. 

The radio source is optically identified with a compact, 
possibly elliptical galaxy located between the two regions 
of high brightness (Richards et al. 1998). It has an iiAB 
magnitude of 24.1, a Gab magnitude of ~26.5, and an H-l-K 
magnitude of 18.7, which corresponds to K~19.0. To date 
no redshift has been measured. However, its faint optical 
magnitude and its red optical to near-infrared colours makes 
it highly unlikely it is at z < 1. Comparison with the K- 
z diagram for powerful radio galaxies (eg. Lilly & Longair 
1984), also makes it most likely that this object is located 
at z > 1. If we assume a lower limit in redshift of z=l, its 
rest-frame 178 MHz power must be higher than PnsMUz ^ 
10^^'^ W Hz~^sr~^, making it at least a few times more 
luminous than HDF123644-I-621133. If it is located at z=1.5. 
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it would have Pitsmhz « 10^^'^ W Hz~^sr~'^. This is such 
a high luminosity that it would be comparable to the most 
powerful FRIs found in the local universe, ie. those near 
the FRI/II division with the optically most luminous host 
galaxies (Ledlow & Owen, 1996). 

Other FRI radio sources in the HDF? 

To identify a radio source at high redshift as one with 
an FRI morphology, the radio observations have to be of 
both sufficient resolution (eg. using MERLIN) and of suf- 
ficient sensitivity to detect the low level extended emis- 
sion. The edge-darkened emission, characteristic for an FRI 
radio source, is in the case of HDF123644-I-621133 and 
HFF123725+621128 only detected a few a level, while these 
sources are two of the brightest objects in the HDF and 
HFFs. Muxlow et al. (1999) have identified 14, much fainter 
sources which appear to have compact two sided emission. 
For a significant fraction of those, fainter emission on larger 
extended scales may have been missed, hampering their clas- 
sification as an FRI. We therefore believe that using the ex- 
isting VLA and MERLIN data, only sources down to mjy 
levels can reliably be classified as FRIs. 



3 THE COSMOLOGICAL EVOLUTION OF 
FRI RADIO SOURCES 

The local space density of FR I radio sources 

We first want to assess the question of how many lumi- 
nous FRI radio galaxies one would expect in an area of 
sky as large as the HDF -|- HFF for a non-evolving cos- 
mological evolution scenario. As a luminosity cut-off we 
take the luminosity of HDF123644-I-621133, which is the 
weaker of the two, independent of the exact redshift of 
HFF123725-^621128. 

Firstly, the local space density of FR I radio sources has 
to be determined. We used the complete subsample of 3CR 
as defined by Laing et al. (1983) for this purpose. It iden- 
tifies 30 FRI radio galaxies in an area of sky with 5 > 10° 
and |6| > 10° (4.24 steradians). Using the 3CR fiux den- 
sity cut-off of 10.9 Jy at 178 MHz and an average spectral 
index of a = —0.75, a radio galaxy with Pitsmhz = lO^'*'* 
W Hz~^sr~^ can be seen out to z=0.042. Four FRI galax- 
ies are found with Pitsmhz > 10^*'* W Hz~^sr~^ in the 
comoving volume of 2.2 x lO^Mpc^ [z < 0.042), which cor- 
responds to a local space density of 170 FRIs per Gpc'^. Note 
that we obtain a very similar value of 200 FRIs per Gpc'^ 
brighter than Pitsmhz > 10^^'* W Hz~^sr~^, using the lo- 
cal luminosity function of Dunlop & Peacock (1990; DP90) 
of steep spectrum radio sources, assuming that all sources 
with Pa.TGHz > lO'* '' W Hz-^sr-i (= Pitsmhz > 25.3 W 
Hz~^sr~^) are FRII, and all the fainter sources are FRIs. 
We will use the latter value in our calculations. 

Implication for a non-evolution scenario for FRI 
radio galaxies 

Hence, what is the chance of finding 2 FR I radio sources in 
an area as large as the HDF and HFF for a non-evolving pop- 
ulation? The area of sky surveyed by the VLA and MERLIN 



covers 10' x 10'. The main parameter involved is the max- 
imum redshift, Zmax, at which the radio sources still would 
have been detected and classified as FRIs. To estimate this 
is a difficult task, not the least due to the uncertain redshift 
of one of the objects. We believe that the limiting factor 
for Zmax for both sources is the low surface brightness level 
of the extended emission in their lobes, which is propor- 
tional to (1 + z)"'*, assuming a spectral index in the lobes 
of a = -1.0. 

In both objects, the extended emission is detected at 
about a 10 o- level, implying that HDF123644-I-621133 could 
have been detected out to z«1.7 at > 3(t. We do not know 
the redshift of HFF123725+621128, but as we argued above, 
we believe it will probably be in the range l<z<1.5~2. If 
it were at z=l, it would have a Zmax of ~1.7. If it were at 
z=1.5 or z—2.0, it would have a Zmax of ~2.4 or ~3.0. We 
therefore do the calculation for Zmax ~ 1.5, 2.0, 2.5 and 3.0. 
The data are given in table with the «max in column 1, 
the comoving volume densities in column 2, the expected 
number of FRIs in the HDF & HFF with Pitsmhz > 24.8 
W Hz~^sr~^ in column 3. The chance, P2, to find two or 
more objects in this volume, if the expected number is p is, 

P>2=e-^{e^~l-p) (1) 

which is derived using a Poisson distribution. This is given 
in column 4 of table It indicates that the presence of two 
FRI radio galaxies in the HDF and HFFs is inconsistent 
(< 1%) with a no-evolution scenario for FRIs, and that 
their comoving space density at z — 1 — 2 was about an 
order of magnitude higher than at present time. 

Constraints on other cosmological evolution models 

Now that we have determined that it is unlikely that FRI 
radio sources undergo no cosmological evolution from z=0 
to z=l, we want to investigate whether FRI and FRII ra- 
dio sources have to be treated as to distinct populations of 
objects. In this section we will compare the data with the 
evolutionary models of DP90. DP90 do not treat the FRIs 
and FRIIs separately; instead they adopt a luminosity- 
dependent cosmological evolution for the total population of 
steep spectrum sources, derived from source counts and red- 
shift information. We will consider two of their models: the 
pure luminosity evolution model and the luminosity /density 
evolution model. 

In both cases, the low redshift luminosity function (LF) 
is fitted using two powerlaw slopes, with the LF being flat- 
ter below the break luminosity and steeper above. In the 
pure luminosity evolution model, the overall shape of this 
LF does not change with cosmological epoch, only the nor- 
malisation in luminosity. Since the slope of the LF is flatter 
at low luminosities, it results in less cosmological evolution 
for the weaker objects. Since we only have to consider ob- 
jects below the the break luminosity, where the slope of the 
luminosity function is a — —0.69, at a specific redshift the 
density is enhanced by a factor lo"'®^ * /. The parameter / 
is the amount of luminosity evolution, which is determined 
by DP90 to be / = 1.26^ - 0.26^^. For the redshift ranges 
involved, the local space density of FRIs would on aver- 
age be enhanced by a factor of 5 — 8, which corresponds to 
P>2 = 2 — 26%. The data are given in table ^ with column 5 
giving the mean density enhancement, column 6 giving the 
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Table 2. Chance of finding 2 FRIs in the HDF and HFF with PirgMHz > 24.8 W Hz '^sr ^ for no evolution, and the pure luminosity 
and luminosity /density evolution models of Dunlop &; Peacock (1990), using difl'erent redshift cut-offs. 
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Figure 2. Average V/Vmax values for FRIs (squares) and FRIIs 
(diamonds), for bins in Alogj^Q(Pi5iiy[jjz) = 0.5. These values are 
taken from Jackson and Wall (1999). Note that no distinction can 
be made between FRIs and FRIIs. 



expected number of FRIs in the HDF & HFF, and column 
7 giving the chance of finding two or more FRIs in this area, 
for the pure luminosity evolution model of DP90. 

The luminosity /density evolution model of DP90 is 
based on their pure luminosity evolution model, but modi- 
fied to allow also the density normalisation to vary with z. In 
a similar fashion as above, but with slightly different param- 
eters, this model would result in an average space density 
enhancement of a factor 4 to 7, which implies a chance of 
finding two or more FRIs in this area of 2—21%, for the 
luminosity/density evolution model of DP90. This is given 
in columns 8—10 of table |^. 

Therefore, in contrast to the non-evolution models, the 
DP90 models, which do not treat FR I and FR II as distinct 
populations, predict surface densities of FRI radio sources 
which are only slightly lower than found in the HDF and 
HFF. 



4 DISCUSSION AND CONCLUSIONS 
A single population scheme? 

The deep optical/radio observations of the HDF and HFF 
have allowed for the first time an estimate of the space den- 
sity of FRI radio sources at high redshifts. The two FRI 
radio sources present in this area of sky indicate that it is 
unlikely that FRI radio sources undergo no cosmological 
evolution. We believe that this may undermine the reason 
for having a "dual population" scheme, in which the FRI 
and FRII radio sources are treated as separate classes of 
objects, with little or no evolution and strong cosmological 
evolution respectively (eg. Jackson & Wall 1999). We want 
to advocate a "single population" scheme, in which the cos- 
mological evolution is a function of radio power, but not de- 
pendent on FR class. Since powerful radio sources undergo a 
stronger evolution than the less powerful ones, and FR Is are 
mainly radio sources of low power, FRIs undergo less evo- 
lution than FRIIs. However, the populations of FRI and 
FRIIs of similar radio power undergo similar cosmological 
evolution. Indeed, the evolution models of Dunlop & Pea- 
cock (1990), which make no distinction between FR classes, 
do fit number-count statistics and redshift distributions well. 
In addition they predict number densities which are compa- 
rable to that found in the HDF and HFF. 

Basically, it is the coincidence of two observations which 
have led to a dual population scheme: 1) the V/Vmax test 
showing that high power radio sources undergo strong space 
density enhancements and low power sources do not, and 
2) the morphological dichotomy between sources of high 
and low radio power. However, the V/Vmax tests for FRI 
and FRII are sensitive over very different redshift ranges. 
Jackson and Wall (1999) show in their figures 2 and 3 the 
mean V/Vmax values as a function of radio power for FRIs 
and FRIIs respectively, indicating an increase in V/Vmax as 
function of radio power. We reproduce these values in figure 
^ using different symbols for the two FR classes. As can be 
seen, no distinction can be made between the mean V/Vmax 
values of FRI and FRII sources at a particular radio power. 
It is due to the fact that FRIs are of lower radio power, and 
therefore only found at the lowest redshifts, that they show a 
lower average V/Vmax- We therefore believe that the results 
of the V/Vmax test actually do not make a dual population 
scheme necessary, and perhaps even argue against one. 

Willott et al. (2001) also introduced a dual population 
scheme in their analysis of complete flux density limited 
samples of 3C, 6C and 7C, which provided an unprecedented 
coverage of the redshift-luminosity plane. Their division be- 
tween low power radio sources with weak emission lines and 
high power radio sources with strong emission lines (not by 
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FR class), results in a space density enhancement of about 
an order of magnitude for the weak objects out to z=l, and 
about three orders of magnitude for the powerful objects 
out to z=2. Using the current data it is difficult to make a 
distinction between this particular scheme and a single pop- 
ulation scheme as proposed here: although the two schemes 
have conceptually different viewpoints, they result quanti- 
tatively in a very similar cosmological evolution for the FR I 
galaxies. Note however, that the evolution as proposed by 
Willott et al. (2001) results in a peculiar luminosity function 
at high redshift, with a prominent 'hump' at the location of 
the break- luminosity at low redshift. 

Treating the population of powerful radio-loud AGN as 
a single class of object would have many benefits, since in 
this way the two FR classes are closely linked. A popular 
paradigm is that the population of radio-loud AGN come 
with a range of jet outputs, of which the more powerful may 
be strong enough to maintain their integrity until they im- 
pact on the intergalactic medium (IGM) in a shock. This 
results in an FRII. However, if the jet is too weak, it may 
dissipate its energy by entraining IGM material, resulting 
in a more turbulent FRI. This may also explain the depen- 
dence of the FRI/FRII radio power divide on the luminosity 
of the host galaxy: in more luminous galaxies, which reside 
in denser environments or which have denser ISM, only jets 
of higher power can keep their integrity. 

If the dual population scheme were to be correct, and 
FRIs and FRIIs undergo different cosmological evolutions, 
it would imply fundamental differences for these two classes 
of object, like the properties of their central engines such 
as black-hole spin or jet composition. As has been argued 
by Gopal-Krishna & Wiita (2000), the existence of radio 
sources, in which the two lobes exhibit clearly different FR 
morphologies, are difficult to reconcile with such models, but 
support explanations for the FR dichotomy based upon jet 
interaction with the medium external to the central engine. 
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